set more off

use dataset, clear

local Z pcteligible2013

keep if abs(x) < 100

su `Z', det
gen Z = `Z' >=r(p50)
replace Z = . if mi(`Z')

su demvoteshare, det
local mymed = r(p50)

preserve

	keep if demvoteshare < `mymed'

	#delimit

	gr tw 
		(sc dpctui20142013 x if Z == 0, msize(tiny) col(gs13))
		(lpolyci dpctui20142013 x if x < 0 & Z == 0, clcol(black) lwid(medthick) ciplot(rline) alpat(solid) alcol(black) alwid(thin))
		(lpolyci dpctui20142013 x if x > 0 & Z == 0, clcol(black) lwid(medthick) ciplot(rline) alpat(solid) alcol(black) alwid(thin))
		,
			ylab(-40(10)10, angle(horiz) nogrid)
			plotregion(style(none))
			subtitle(, fcol(none) lwid(none))
			xline(0)
			ylab(, angle(horiz))
			xtitle("Distance from nearest state border")
			legend(off)
			ytitle("Change in pct. uninsured")
			title("Low eligibility")
			legend(off)
			nodraw
			name(g1, replace)
			;		
		
	#delimit cr

	#delimit

	gr tw 
		(sc dpctui20142013 x if Z == 1, msize(tiny) col(gs13))
		(lpolyci dpctui20142013 x if x < 0 & Z == 1, clcol(black) lwid(medthick) ciplot(rline) alpat(solid) alcol(black) alwid(thin))
		(lpolyci dpctui20142013 x if x > 0 & Z == 1, clcol(black) lwid(medthick) ciplot(rline) alpat(solid) alcol(black) alwid(thin))
		,
			ylab(-40(10)10, angle(horiz) nogrid)
			plotregion(style(none))
			subtitle(, fcol(none) lwid(none))
			xline(0)
			ylab(, angle(horiz))
			xtitle("Distance from nearest state border")
			legend(off)
			ytitle("Change in pct. uninsured")
			title("High eligibility")
			legend(off)
			nodraw
			name(g2, replace)
			;
		
	#delimit cr


	#delimit

	gr tw 
		(sc dpctui20152013 x if Z == 0, msize(tiny) col(gs13))
		(lpolyci dpctui20152013 x if x < 0 & Z == 0, clcol(black) lwid(medthick) ciplot(rline) alpat(solid) alcol(black) alwid(thin))
		(lpolyci dpctui20152013 x if x > 0 & Z == 0, clcol(black) lwid(medthick) ciplot(rline) alpat(solid) alcol(black) alwid(thin))
		,
			ylab(-40(10)10, angle(horiz) nogrid)
			plotregion(style(none))
			subtitle(, fcol(none) lwid(none))
			xline(0)
			ylab(, angle(horiz))
			xtitle("Distance from nearest state border")
			legend(off)
			ytitle("Change in pct. uninsured")
			title("Low eligibility")
			legend(off)
			nodraw
			name(g3, replace)
			;
		
	#delimit cr

	#delimit

	gr tw 
		(sc dpctui20152013 x if Z == 1, msize(tiny) col(gs13))
		(lpolyci dpctui20152013 x if x < 0 & Z == 1, clcol(black) lwid(medthick) ciplot(rline) alpat(solid) alcol(black) alwid(thin))
		(lpolyci dpctui20152013 x if x > 0 & Z == 1, clcol(black) lwid(medthick) ciplot(rline) alpat(solid) alcol(black) alwid(thin))
		,
			ylab(-40(10)10, angle(horiz) nogrid)
			plotregion(style(none))
			subtitle(, fcol(none) lwid(none))
			xline(0)
			ylab(, angle(horiz))
			xtitle("Distance from nearest state border")
			legend(off)
			ytitle("Change in pct. uninsured")
			title("High eligibility")
			legend(off)
			nodraw
			name(g4, replace)
			;
		
	#delimit cr

	gr combine g1 g2, rows(1) title("Change in uninsured 2014-2013") name(G1, replace) 
	gr combine g3 g4, rows(1) title("Change in uninsured 2015-2013") name(G2, replace) 
	gr combine G1 G2, rows(2) name(GA, replace) title("Republican counties")
	gr export "_output/figureA14.pdf", replace
	
restore

preserve

	keep if demvoteshare > `mymed'

	#delimit

	gr tw 
		(sc dpctui20142013 x if Z == 0, msize(tiny) col(gs13))
		(lpolyci dpctui20142013 x if x < 0 & Z == 0, clcol(black) lwid(medthick) ciplot(rline) alpat(solid) alcol(black) alwid(thin))
		(lpolyci dpctui20142013 x if x > 0 & Z == 0, clcol(black) lwid(medthick) ciplot(rline) alpat(solid) alcol(black) alwid(thin))
		,
			ylab(-40(10)10, angle(horiz) nogrid)
			plotregion(style(none))
			subtitle(, fcol(none) lwid(none))
			xline(0)
			ylab(, angle(horiz))
			xtitle("Distance from nearest state border")
			legend(off)
			ytitle("Change in pct. uninsured")
			title("Low eligibility")
			legend(off)
			nodraw
			name(g1, replace)
			;		
		
	#delimit cr

	#delimit

	gr tw 
		(sc dpctui20142013 x if Z == 1, msize(tiny) col(gs13))
		(lpolyci dpctui20142013 x if x < 0 & Z == 1, clcol(black) lwid(medthick) ciplot(rline) alpat(solid) alcol(black) alwid(thin))
		(lpolyci dpctui20142013 x if x > 0 & Z == 1, clcol(black) lwid(medthick) ciplot(rline) alpat(solid) alcol(black) alwid(thin))
		,
			ylab(-40(10)10, angle(horiz) nogrid)
			plotregion(style(none))
			subtitle(, fcol(none) lwid(none))
			xline(0)
			ylab(, angle(horiz))
			xtitle("Distance from nearest state border")
			legend(off)
			ytitle("Change in pct. uninsured")
			title("High eligibility")
			legend(off)
			nodraw
			name(g2, replace)
			;
		
	#delimit cr

	#delimit

	gr tw 
		(sc dpctui20152013 x if Z == 0, msize(tiny) col(gs13))
		(lpolyci dpctui20152013 x if x < 0 & Z == 0, clcol(black) lwid(medthick) ciplot(rline) alpat(solid) alcol(black) alwid(thin))
		(lpolyci dpctui20152013 x if x > 0 & Z == 0, clcol(black) lwid(medthick) ciplot(rline) alpat(solid) alcol(black) alwid(thin))
		,
			ylab(-40(10)10, angle(horiz) nogrid)
			plotregion(style(none))
			subtitle(, fcol(none) lwid(none))
			xline(0)
			ylab(, angle(horiz))
			xtitle("Distance from nearest state border")
			legend(off)
			ytitle("Change in pct. uninsured")
			title("Low eligibility")
			legend(off)
			nodraw
			name(g3, replace)
			;
		
	#delimit cr

	#delimit

	gr tw 
		(sc dpctui20152013 x if Z == 1, msize(tiny) col(gs13))
		(lpolyci dpctui20152013 x if x < 0 & Z == 1, clcol(black) lwid(medthick) ciplot(rline) alpat(solid) alcol(black) alwid(thin))
		(lpolyci dpctui20152013 x if x > 0 & Z == 1, clcol(black) lwid(medthick) ciplot(rline) alpat(solid) alcol(black) alwid(thin))
		,
			ylab(-40(10)10, angle(horiz) nogrid)
			plotregion(style(none))
			subtitle(, fcol(none) lwid(none))
			xline(0)
			ylab(, angle(horiz))
			xtitle("Distance from nearest state border")
			legend(off)
			ytitle("Change in pct. uninsured")
			title("High eligibility")
			legend(off)
			nodraw
			name(g4, replace)
			;
		
	#delimit cr

	gr combine g1 g2, rows(1) title("Change in uninsured 2014-2013") name(G1, replace) 
	gr combine g3 g4, rows(1) title("Change in uninsured 2015-2013") name(G2, replace) 
	gr combine G1 G2, rows(2) name(GB, replace) title("Democratic counties")
	
	
	gr export "_output/figureA15.pdf", replace

restore
gr close
